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T^ 1 INTRODUCTION 

<D . 

^ ■ In this paper, I intend to motivate supersymmetric grand unified theories (SUSY 

GUTs), briefly explain an extension of the standard model based on them and present 
a calculation performed using certain properties of some SUSY GUTs to constrain the 
H ' available parameter space. 

Why GUTs? 

Much work has been done on the running of the gauge couplings in the standard 
model, as prescribed by the renormalisation group. Amazingly, when the couplings cui, 
a?2 and a.3 were run up to fantastically high energies ~ 10 14 GeV, they seemed to be 
converging 1 to one value. This is a feature naturally explained by many GUTs such 
as SU(5) 2 ' 3 and reflects the fact that the strong, weak and electromagnetic forces seen 
today are different parts of the same grand unified force. It was realised that GUTs 
could also provide relations between the masses of the observed fermions, the structure 
and hierarchy of which are as yet unexplained. Despite these attractive features, several 
problems arose which detracted from the idea. 

Unfortunately, the three couplings do not quite converge by ~ 7a, and many GUTs, 
notably SU(5), predict proton decay much faster than the lower experimental bounds. 
Also incredible fine tuning is required for the so-called 'hierarchy problem'. This stems 
from the fact that M\y changes through radiative corrections of order the new physics 
scale (Fig.l) , say the Planck mass ~ 10 19 GeV, if there is no new physics at smaller 
energies. My/ is therefore unstable to the corrections and vast cancellations in the 
couplings are required to motivate the correct phenomenology. 




Figure 1: One loop corrections to m\. The first diagram gives a ~ Mp ; contribution. 




Figure 2: One loop triple hypercharge boson anomally. 



Why SUSY? 

Supersymmetry is an extra symmetry relating fermions and bosons, and so provides 
some explanation of how particles of differing spin should be related to one another. 
In an unbroken supersymmetric theory, each fermion has a degenerate bosonic partner. 
Of course, these so called superpartners are not observed, so that if supersymmetry was 
ever the correct theory, it must have been broken. However, with the introduction of 
superpartners at some rough energy scale Msusy, the renormalisation group running 
of the gauge couplings changes. The coupling constants are now seen to meet at a scale 
~ 10 16 GeV, as reflected by the correct sin 2 9 W prediction 4 . My/ becomes stabilised 
because supersymmetry induces cancellations between the bosonic and fermionic loop 
corrections to the mass. The quadratic divergences induced by the loop corrections now 
add to zero and one is left with merely logarithmic divergences. 

THE MINIMAL SUPERSYMMETRIC STANDARD MODEL (MSSM) 

The MSSM is a minimal extension of the standard model into supersymmetry. In the 
model, every particle of the standard model has a superpartner associated with it that 
transforms identically under the standard model gauge group but have spin different 
by |. So for example, each quark has a scalar "squark" superpartner, the gluons have 
"gluinos" etc. At first sight however, the model has a U(l) y gauge anomaly. This 
originates from the diagram with three B gauge bosons connected to an internal loop 
through which any fermions may run (cf Fig.2) and the counter term to it would destroy 
gauge invariance. The diagram is proportional to J2i(Xi/2) 3 where i runs over all active 
fermions. Through the hypercharge assignments, this cancels in the standard model 
but in the MSSM the superpartner of the Higgs called the Higgsino with Y = 1 may 
run around the loop. To cancel this effect, a second Higgs H 2 must be introduced which 
transforms in the same way to Hi except for having Y — — 1. 



The new Higgs must also develop a vev v 2 to give masses to up quarks and the two 
vevs are related by 

tan/3 = — (1) 
vi 

where v\ + v\ — v 2 and v = 246 GeV, the measured vev of the standard model. 
In chiral superfield form, the superpotential looks like 

Wmssm = UQH 2 u c + DQH x d c + ELH ie c + ^H 1 H 2 (2) 

where U, D and E are the up, down and charged lepton Yukawa matrices respectively 
and all gauge and family indices have been suppressed. 

One possible problem with this superpotential is the dimensionful parameter fi. \x 
needs to be ~ Mz to give the right electroweak symmetry breaking behaviour whereas 
one would expect it to be of order of the new physics scale M GUT . One solution to 
this problem is described in the Next to Minimal Supersymmetric Standard Model 
(NMSSM). 

The NMSSM 

The fi term in Eq.|2] is replaced by XNH1H2 where N is a gauge singlet and therefore 
doesn't affect the coupling constant unification. In certain supergravity models, N 
develops a vev naturally of order Mz and so the /x term is generated without having 
to put fj, in "by hand." The superpotential now has a discrete Peccei-Quinn symmetry 
which leads to phenomenologically unacceptable low energy axions and so a term — |iV 3 
is added which breaks it. [] 



GUTS WITH YUKAWA UNIFICATION 

GUTs can quite naturally provide Yukawa unification relations between the quarks 
and/or leptons. For example in SU(5), the right handed down quarks and conjugated 
lepton doublet lie in a 5 representation. When a mass term ~ 5*5j is formed, the 
Yukawa relation 

K(M GUT ) = X T (M GUT ) (3) 

applies. Also in SO(10), the whole of one family and a right handed neutrino is con- 
tained in one 16 representation, leading to triple Yukawa unification, where the top, 
bottom and charged lepton Yukawa couplings are equal at the GUT scale. 

These relations can be used to constrain the parameter space of m t and tan j3, which 
has been done for the MSSM 5 . Our idea was to repeat this calculation for the NMSSM, 
to see how much the viable parameter space changes in the model. 



THE CALCULATION 

The basic idea is to choose some tan /3 and m t and run and A r up to M GUT ~ 10 16 
GeV. Then, to some arbitrary accuracy, one can determine whether the GUT relation 
EqfJ holds. If it does, then SU(5) and the other Yukawa unifying extensions of the 
standard model are possible on this point in parameter space. The procedure is iterated 
over all reasonable values of tan/? and m t . The calculation is presented in more detail 
in Ref.6. 



A and k are merely coupling constants 
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Figure 3: Viable Range of Parameter Space For as(Mz) 
values are quoted at m t . 



0.11, m b = 4.25 GeV. A and k 



Starting Point M z - 

We use the definitions of the gauge couplings at M z : a 1 " 1 (M 2 ) = 58.89, a^ l (M z ) = 
29.75 and a^ l (M z ) = 0.11 ±0.01. The first two gauge couplings are determined 
accurately enough for our purposes whereas the third needs to be used as a parameter, 
on account of its large uncertainty. 

In order to convert masses of quarks to Yukawa couplings, we simply need to read 
them off the potential Eq.0 at some energy scale (taken here to be mt): 



where 



A 4 (m t ) 
At (m t ) 
X T (m t ) 
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v sin (3 
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(4) 
(5) 
(6) 

(7) 



rrif (mt) ' 

Note that whereas the mt referred to here is always the running one, it can be related 
to the physical mass by 5 



m phys _ m ^ 1 + + q 



(8) 



To determine 7] b and r] T , the masses are run up from the on shell mass to m t using 
effective 3 loop QCD ® 1 loop QED 7 > 8 ' 9 > 10 . Note that these factors will depend of 
m b = 4.25 ± 0.15 GeV and a 3 (M z ). m t is assumed to be the rough energy scale 
when the whole supersymmetric spectrum kicks in. While being unrealistic, trials 
with Msusy — 1 TeV show only a few percent deviation from the predictions with 
Msusy = rrit- So, having determined the gauge and relevant Yukawa couplings at m t , 




Figure 4: Viable Range of Parameter Space For as(Mz) 
m b = 4.1-4.4 GeV. The left most lines are for X b = 0.9A r . 



0.12 and experimental bounds of 



we need RG equations to run them up to Mqut in the NMSSM. To derive these, we 
used results from a general superpotential 11 to obtain 



16^ 



16n 



16tt z 



16tH 



dt 
,d\b 
dt 
>dX T 

~dT 

,dX 



IQti 2 



dt 
,dk 

dt 



A, 
h 

X 
6A- 



6A 2 + A 2 + A 2 - ( f/, + 3gj + fgl 



6A 2 + A 2 + A 2 + A 



7_ 

15 



9 2 i + ^1 + y^ 3 2 ) 



A 2 + 3A^ + A 2 



gl + Zgl 



4A 2 + 2k 2 + 3X 2 T + 3A^ + 3A 
A 2 + A; 2 



:^i 2 + 3^ 2 



(9) 



in the limit that the lighter two families have negligible contributions (a very good 
approximation) . 

The Yukawa couplings can now be run from m t to 10 16 GeV using numerical tech- 
niques. The parameters A and k particular to the NMSSM are unconstrained at m t so 
they are merely varied for different curves. 

Our results are displayed in Fig. as contours in the tan f3 — m t plane consistent 
with Eq.|. We take a 3 (M z ) = 0.11, m b = 4.25GW and the NMSSM parameters 
A (rat) and fc(rat) as indicated. The MSSM contour is shown for comparison and is 
indistinguishable from the NMSSM contour with A (ra t ) = 0.1 and k (m t ) = 0.5. In fact 
our plot for the MSSM based on 1-loop RG equations is very similar to the 2-loop result 
in ref.5. The deviation of the NMSSM contours from the MSSM contour depends most 
sensitively on A (rat) rather than k(mt). Two of the contours are shortened due to either 
A or k blowing up at the GUT scale. For A (rat) = 0.5, k (m t ) = 0.5, no points in the 
rat — tan/5 plane are consistent with Eq.||| Yukawa unification, while for A (rat) = 0.1, 
k (mt) = 0.1 — 0.5 the contours are virtually indistinguishable from the MSSM contour. 
In general we find that for any of the current experimental limits on a 3 and ra&, the 
maximum value of A(rat) or k(mt) is ~ 0.7 for a perturbative solution to Eq||. 

Fig.|] shows the effects of particle thresholds, which can modify Eq|| to Xb = 0.9A T . 
Our treatment does not treat supersymmetric or heavy thresholds exactly and so some 



sort of corrections like those shown are expected. The curves are at ctg(Mz) = 0.12 
and nib = 4.1-4.4 GeV to illustrate that uncertainties in these quantities make a large 
difference to the parameter space. These uncertainties are much bigger than those 
associated with the NMSSM, and so the MSSM and NMSSM would be practically 
indistinguishable given the parameters m t and tan f3. 

Other Yukawa Parameters 

The next useful step is to notice that Eqs.pl are all of the form 
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where M" 1 and c" are constants supplied by the relevant RG equation. When the (3 
function 

dg% = hgf_ (u) 

dt 16tt 2 1 ' 

is inserted, and the RG equations are reparameterised in terms of the flow and not the 
trajectory of the solutions, we obtain 



A (M; 
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K(M GUT ) 

where 



Cexp(-]TM^), (12) 



C a - TT ( a{McuT) \ 2 " (13) 

^ -jiWMsusy)) (13) 



contains all the information about the gauge couplings and 

An(M GUT ) 



167T 2 Jln(MsusY) 



Xjdt (14) 



concerns the Yukawa couplings. 

With this formulation, the running of the physically relevant Yukawa eigenvalues 
and mixing angles can be expressed in simple terms as shown below, 
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(15) 



with identical scaling behaviour to V c b of V u b, Vt s , Vtd- To a consistent level of ap- 
proximation V us , V ud , V cs , V cd , V tb , A„/A c , Xd/X s and Ag/A^ are RG invariant. The CP 
violating quantity J scales as V c \. Eqs. [14| also apply to the NMSSM since the 
extra A and k parameters cancel out of the RG equations in a similar way to the gauge 
contributions as can easily be seen from Eq||. The only difference to these physically 
relevant quantities is therefore contained in I T , lb and I t . 
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Figure 5: h integrals for as(Mz) = 0.11, m& = 4.25 GeV. 

These Jj integrals are shown in Fig.|5| and the NMSSM results are the upper lines of 
each pair, and it is clear that the deviation between the two models is small again. 

We emphasise that the results of the Ii integrals shown in Fig.|5| play a key role 
in determining the entire fermion mass spectrum via the scaling relations of Eq.|I~5|. 
The small deviation between the NMSSM and the MSSM results compared to the 
experimental uncertainties means that the recent GUT scale texture analyses of the 
quark mass matrices which were performed for the MSSM are equally applicable to the 
NMSSM. For example, the recent Ramond, Roberts and Ross (RRR) 12 texture analysis 
is also based upon Eq|| and assumes a Georgi-Jarlskog (GJ) 13 ' 14 ansatze for the charged 
lepton Yukawa matrices, although their results in the quark sector are insensitive to 
the lepton sector. It is clear that all the RRR results are immediately applicable 
to the NMSSM since the only difference between the two models enters through the 
scaling integrals Ij whose deviation we have shown to be negligible compared to the 
experimental errors. 

CONCLUSIONS 

We have discussed the unification of the bottom quark and tau lepton Yukawa 
couplings within the framework of the NMSSM. By comparing the allowed regions of 
the m t -tan/3 plane to those in the MSSM we find that over much of the parameter 
space the deviation between the predictions of the two models which is controlled by 
the parameter A is small, and always much less than the effect of current theoretical and 
experimental uncertainties in the bottom quark mass and the strong coupling constant. 
We have also discussed the scaling of the light fermion masses and mixing angles, and 
shown that to within current uncertainties, the results of recent quark texture analyses 12 
performed for the minimal model also apply to the next-to-minimal model. There are 
however two distinguishing features of the NMSSM. Firstly, the scaling of the charged 
lepton masses will be somewhat different, depending on A and k. Although this will not 
affect the quark texture analysis of RRR, it may affect the success of the GJ ansatze 13 ' 14 
for example. Secondly, the larger tan/3 regions may not be accessible in the NMSSM 
for large values of A and k, so that full Yukawa unification may not be possible in this 
case. 
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